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SUMMARY 


\ 


The  aerofoil  NLR  7301,  16.5%  thick,  was  designed  at  NLR  by  their  generalised 
hodograph  method  to  give  supercritical  (inviscid)  shock-free  flow  at  the  design 
condition  M^  = 0.72,  CL  = 0.6.  The  paper  describes  a theoretical  study  of  this 
aerofoil  by  the  latest  form  of  the  RAE  ' VGK'  method  for  calculating  viscous 
effects  in  two-dimensional  transonic  flows.  First,  comparisons  are  given  with 
experiments  in  the  NLR  Pilot  Tunnel  at  a low  Reynolds  number  (2  x 10  ).  Large 
discrepancies  are  shown  near  the  design  condition  and  at  higher  Mach  numbers, 
which  are  ascribed  to  a laminar  (transitional)  shock  wave/boundary  layer  inter- 
action in  the  transition-free  experiments;  this  causes  an  apparent  weakening  of 
the  shock  wave  and  a spuriously  low  level  of  the  measured  drag.  Next,  calculations 
were  made  at  two  higher  Reynolds  numbers,  with  transition  fixed  near  the  leading 
edge:  10^*  to  simulate  conditions  in  a medium-sized  wind  tunnel;  and  5 x 10  , 
typical  of  full-scale  conditions.  Favourable  scale  effects  are  predicted  on  drag- 
rise  Mach  numbers,  for  Cl  * 0.5,  and  on  separation  onset  at  all  lift  coefficients. 
It  is  concluded  that  the  aerofoil  should  have  an  excellent  performance  at  high 
Reynolds  numbers.  _ v , i 
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1 INTRODUCTION 

The  aerofoil  NLK  7301  is  a lb. 5%  thick  supercritical  aerofoil,  designed  by 

the  generalised  hodograph  method  of  Boerstoel  and  Huizing*;  it  was  selected  by 

Nl.R  for  inclusion  in  the  'data  base'  of  experimental  information  supplied  to  the 

AGARI)  FOP  WG4 . The  present  Report  describes  an  extensive  theoretical  study  of 
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this  aerofoil  by  the  method  of  Collyer~  and  hock  ’ ; comparisons  with  experiments 
in  the  NLR  Pilot  Tunnel  (at  a Reynolds  number  of  2 * 10  ) are  included,  but  most 
of  the  calculations  have  been  made  at  higher  Reynolds  numbers  - 1 0 ^ to  represent 
conditions  in  a larger  wind  tunnel,  and  5 * 1 0 ^ to  represent  lull  scale  conditions. 

2 THK  NUNKR1CAL  MKT HOD 

The  method  of  calculation,  described  briefly  in  Ref  4 (p  4t>9)  and  fully 
in  Ref  2,  is  an  adaption  of  the  inviscid  finite  difference  method  of  Carabedian 
and  Korn"*  (G  & K),  combining  it  with  the  'lag-entrainment'  method  of  Green  of  <j2h 
for  calculating  the  boundary  layer  and  wake  to  allow  fully  for  viscous  effects; 
the  'transpiration'  model  of  the  displacement  effect  is  used  to  provide  an  inner 
boundary  condition  for  the  inviscid  calculations,  both  on  the  aerofoil  surface 
and  in  the  wake.  The  method  has  been  recently  improved  by  incorporating  a variant 
of  Jameson's*'  ’quasi-conservative’  (Q-G)  difference  scheme;  a parameter  \ has 
been  introduced  in  the  ’shock-point'  difference  operator  which  is  such  that  A » 0 
gives  a non-conservative  (N-C)  scheme,  A « I gives  the  Q-C  scheme,  while  inter- 
mediate values  naturally  lead  to  intermediate  results.  Since  it  has  been 
generally  noted  that  neither  the  'N-G*  or  'Q-C'  extremes  give  uniformly  good 
agreement  with  experiment  (when  viscous  effects  are  allowed  for),  it  is  not 
surprising  that  a value  of  \ of  about  0.5  does  generally  give  better  agreement 
with  experiment:  see  Ref  1 for  example.  This  value  of  A has  been  used  in  all 
the  calculations  described  below.  Otherwise  the  standard  parameters  of  the 
method  have  been  used*:  results  are  all  for  the  'fine'  (IbO  * 30)  grid,  making 
normally  50  iterations  on  the  'coarse'  (80  * 15)  grid  followed  by  about  200  on 
the  fine  grid. 

1 GALUULATIONS  FOR  1NV1SGID  FLOW 

The  aerofoil  is  shown  in  Fig  I.  Calculations  of  the  pressure  distribution 
tor  the  design  condition  (M  » 0.721,  a = - 0.2°,  G,  = 0.599)  are  shown  in  Fig  2; 
the  lull  line  being  the  'exact'  (shock-free)  solution  obtained  by  the  design 
method  and  the  dashed  line  that  obtained  by  the  present  direct  method.  Agreement 
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* The  'artificial  viscosity’  parameter  < (IT)  is  set  at  0.8  throughout. 
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is  excellent  except  near  the  rear  sonic  point,  and  there  the  shock  predicted 
by  the  direct  method  is  so  weak  that  a drag  coefficient  of  only  0.00016  is 
calculated.  Fig  3a&b  show  how  the  pressure  distribution  develops  at  a * -0.20 
as  the  free  stream  Mach  number  (M  ) is  increased  from  0.68  through  the  design 

oo 

Mach  number  (0.72)  up  to  0.73.  As  expected  a shock  of  moderate  strength  appear 
near  the  leading  edge  (M  - 0.68)  which  moves  back  as  M is  increased,  growing 

OO  do 

weaker  while  a second  (also  weak)  shock  appears  at  — — 0.6  (M^  “ 0.715).  At 
the  design  condition  these  two  coalesce  to  form  one  weak  shock,  which  grows 
rapidly  in  strength  as  is  further  increased.  A rather  similar  pattern  is 

seen  in  Fig  4a&b,  which  show  how  the  presence  distribution  varies  with  o (from 
-0.5°  to  +0.1°)  at  the  design  Mach  number  (0.72).  As  is  usual  with  a super- 
critical aerofoil  of  this  type,  quite  small  changes  in  either  M or  a from 
their  design  values  came  appreciable  changes  in  the  pressure  distribution,  but 
Fig  5 shows  that  for  values  of  a or  below  the  design  values  the  drag 

'creep'  is  small,  not  exceeding  about  0.0005  in  Cp  . Above  the  design  point, 
the  drag  increases  rapidly. 

4 VISCOUS  THEORY:  COMPARISONS  WITH  EXPERIMENT 

A model  of  the  aerofoil,  of  0.18m  chord,  has  been  tested  in  the  NLR  Pilot 
tunnel  (0.55  m x 0.42  m)  at  a Reynolds  number  of  about  2 x 10^.  Because  of  the 

low  Reynolds  number,  the  dilemma  of  whether  to  test  with  transition  'fixed'  or 

'free'  was  more  acute  than  usual;  in  fact  tests  of  both  types  have  been  made, 
though  as  we  shall  see  later  each  type  has  its  particular  defects. 

First,  we  note  from  Fig  6 the  very  large  reduction  in  lift  predicted  by 

allowing  for  viscous  effects  at  the  design  condition  (M  * 0.72,  a = - 0.2°); 

the  lift  coefficient  is  reduced  from  0.6  to  0.3,  almost  exactly  halved.  (Here, 
the  position  of  transition  was  taken  to  be  at  the  suction  peak  on  the  upper 
surface  (x/c  — 0.07),  and  at  x/c  = 0.3  on  the  lower  surface.) 

We  consider  next  the  combination  of  M^  and  a which  appears  at  first 
sight  to  lead  to  an  experimental  distribution  (transition  free)  most  similar  to 
that  for  the  inviscid  design  condition.  This  is  at  = 0.747,  a = 0.85° 
(uncorrected)  and  is  shown  in  Fig  7.  In  an  endeavour  to  obtain  a good  match 
between  viscous  theory  and  experiment,  the  following  procedure  has  been  used  in 
this  and  other  cases.  First  the  presence  ratio  p/H^  over  the  first  20%  chord 
on  the  upper  surface  is  matched  by  varying  a , keeping  M^  fixed.  Next  M(  is 
varied:  as  is  well  known,  the  values  of  p/H^  on  the  upper  surface  near  the 
leading  edge  suction  peak  are  almost  invariant  with  M^  when  the  flow  is  super- 
critical so  that  a match  with  the  lower  surface  pressure  distribution  over  the 
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trout  part  ot  t ho  aerofoil  can  bo  obtained  without  disturbing  that  already 
achieved  on  the  upper  surlaee.  In  tact,  the  value  ot  M arrived  at  in  this 
way  was  tound  to  be  very  close  to  the  nominal  value,  suggesting  that  uo  blockade 
correction  is  rei|uired:  somewhat  surprising  in  view  ot  the  rathet  large  (10") 
open  area  ratio  ot  the  slotted  walls,  which  would  suggest  that  a negative  block- 
age  correct  ion  might  be  needed.  l'ho  position  ot  transition  was  specified  to  be 
as  far  back  as  possible  on  both  surfaces.  On  the  upper  surface  this  is  at  the 
suction  peak  just  before  the  shock  (x/c  “ 0.r>8);  on  the  lowci  surface  it  proved 
impossible  to  get  the  transition  point  further  at  t than  x/c  ■ 0.1/’*,  In  this 
way  very  good  agreement  with  experiment  over  the  trout  third  o!  the  aerofoil  is 
achieved;  but  further  aft  the  situation  is  very  different.  the  most  noticeable 
discrepancy  is  on  the  upper  surface:  here  theory  suggests  that  there  should  be  a 
marked  shock  wave  at  x/c  " O.b,  strong  enough  to  cause  about  -t0  1 counts'  ot  wave 
drag;  while  in  the  experiment  there  is  apparently  no  shock  at  all.  What  has 
almost  certainly  happened  is  that  a laminar  shock-wave  .'boundary- 1 ivet  interaction 
has  taken  place,  which  has  smeared  out  the  pressure  rise  through  the  shock  to  such 
anextont  that  it  has  become  unite  smooth.  This  phenomenon  has  also  caused  the 
wave  d rag  associated  with  the  shock  t o be  negligible:  tin*  value  ot  tneusui ed 

tbv  a pitot  transverse  in  the  wake  0.8  chords  downstream  ot  the  trailing  edge! 
is  only  0.008s  .is  opposed  to  the  calculated  value  0.01.M,  It  is  signit leant 
that  the  purely  'viscous1  part  ot  the  calculated  drag  l neglecting  wave  drag! 
is  0.008-’,  agreeing  closely  with  the  experimental  value.  Making  the  reasonable 
assumption  that  transition  from  laminar  to  turbulent  tlow  is  taking  place  in 
the  experiment  at  about  tin*  position  (x/c  - O.tO  assumed  in  the  ca  l cu l at  ions , wo 
can  deduce  that  little  wave  drag  is  present  in  the  experimental  measurements; 
whereas  the  calculat ions  suggest  that  there  would  be  about  40  ' counts'  ot  wave 
drag  if  the  boundary  layer  were  turbulent  at  its  interaction  with  the  shock  wave. 
On  tin*  lower  surface,  agreement  between  theory  and  experiment  is  better,  though 


* The  reason  for  this  is  that  Tliwaites'  method,  which  is  used  to  calculate  the 
laminar  part  ot  the  boundary  layer,  breaks  down  (indicating  a lamina!  separa- 

t>-  d0  . . 

tionl  when  the  parameter  - V “ - — - exceeds  the  critical  value  0,08.  It 
' v dx 

a transition  point  is  specified  such  that  this  critical  value  is  reached,  then 
the  program  automat ical ly  moves  the  transition  position  forward  until  a non- 
critical  situation  occurs  This  may  easily  happen  when  the  specified  position 
is  far  back  on  the  chord,  even  in  an  apparently  favourable  pressure  gradient, 
because  the  disturbance  to  the  pressure  distribution  at  transition  may  induce 
a local  gradient  sufficient  to  trigger  the  process  (d  being  relatively 

large,  so  that  only  a small  value  ot  dl'/dx  is  needed  to  cause  \ t o teach 
its  or i t i ca I va lue 1 . 
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there  are  signs  of  a laminar  separation  bubble  near  x/c  - 0.5  which  cannot  of 
course  be  predicted  by  the  theory. 

The  results  for  the  same  case,  but  with  transition  fixed  at  x/c  = 0.3  on 
both  surfaces,  are  shown  in  Fig  8.  With  a normal  calculation  (full  line)  having 
these  transition  positions,  the  agreement  with  experiment  on  the  upper  surface 
is  again  poor  behind  x/c  = 0.3,  while  over  the  rear  part  of  the  lower  surface 
the  measured  pressures  are  much  lower  than  those  predicted,  suggesting  that  the 
roughness  band  has  thickened  the  boundary  layer  sufficiently  to  bring  it  near  to 
separation  at  about  707!  chord.  Some  support  to  this  hypothesis  is  provided  by 
some  additional  calculations  made  with  an  increment  in  momentum  thickness 
(A0  ) artificially  added  at  the  transition  point  (on  both  surfaces),  shown  by 
the  dashed  lines;  clearly  the  pressures  are  now  nearer  the  measured  ones  on 
both  surfaces.  Note  that  a double  shock  pattern  occurs  on  the  upper  surface, 
which  is  not  predicted  by  the  theory. 

Figs  9 to  12  compare  theory  with  experiment  (transition  free)  at  a constant 
value  of  a (0.85°)  in  the  latter,  through  the  set  of  Mach  number  values  0.6, 

0.7,  0.72  and  0.77.  At  the  three  lower  Mach  numbers  agreement  is  quite  good, 
except  in  the  region  of  laminar  separation  on  the  lower  surface  mentioned 
previously,  and  also  over  the  last  \0Z  chord  on  the  upper  surface,  where  the 
present  method  is  usually  found  to  overestimate  the  pressure  rise  to  the  trai ling- 
edge  slightly.  However  at  M^  — 0.774  (Fig  12)  there  are  again  signs  of  a marked 

laminar  interaction  in  the  experiment,  lit  was  not  possible  to  reach  this  Mach 

number  in  the  calculations;  even  at  M - 0.765,  boundary  layer  separation  is 
predicted  from  the  foot  of  the  shock  to  the  trailing  edge,  while  at  higher  values 
of  M the  method  diverged.) 

CO 

Figs  1 3 to  15  compare  results  at  constant  M (0.745)  as  «_  varies  from 

O')  00  T 

-2  through  0 to  +2  , At  a^.  “ - 2 we  have  chosen  the  'transition  fixed' 
conditions  in  the  experiment.  On  the  upper  surface  agreement  is  good  but  on  the 
lower  surface  it  is  probable  that  separation  is  imminent  (note  the  high  value  of 

H , 2.2,  predicted  at  x/c  • 0.8),  provoked  by  the  strong  shock  wave  at 

x/c  • 0.45  and  the  subsequent  adverse  pressure  gradient.  At  a = 0 (transition 
free)  agreement  is  fair  on  both  surfaces  (Fig  14);  but  at  = + 2°  (transition 
fixed)  (Fig  15)  it  was  impossible*  to  reach  a sufficiently  high  value  of  a in 


* Because  the  calculations  diverged  (owing  to  excessively  high  values  of  H ). 
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the  calculations  to  match  with  the  experimental  pressures  over  the  first  251 
chord  - even  so  the  agreement  is  moderately  good  on  both  surfaces  and  in  particu- 
lar it  is  clear  that  a laminar  interaction  has  been  avoided  in  the  experiment 
by  fixing  transition  ahead  of  the  shock. 

Our  final  comparison  with  experiment  is  in  Kig  lb,  which  shows  measured  and 

predicted  values  of  for  the  'transition  tree’  condition  with  a.^.  * 0.8 Sl 

(taking  the  position  of  transition  in  the  calculations  as  far  back  on  the  chord  as 

the  method  will  allow*;  as  marked  in  the  appropriate  figures).  1 1 wi 1 1 be  seen  t hat  t he 

theory  slightly  overestimates  the  drag  from  M - O.b  to  0.7  (usually  for  other 

9 . 

aerofoils,  with  transition  fixed  near  the  leading  edge,  it  is  found  to  uttdei'- 
estimate  by  about  the  same  amount).  However,  the  most  marked  discrepancy 

is  that  the  ’bucket'  in  the  experimental  measurements  at  - 0.747  is  replaced 
in  the  theory  by  a rapid  rise  in  drag  at  about  the  same  Mach  number.  This 
strongly  suggests  that  the  laminar  shock  wave/boundary  layer  interaction  present 
in  the  experiment  is  producing  unrealistically  low**  values  of  the  drag,  parti- 
cularly at  the  higher  Mach  numbers  where  the  shock  is  far  back  along  the  chord. 

Note  also  the  dotted  curve  marked  ' C„  . ’ (twice  the  momentum  thickness  cal- 

D viscous 

culated  at  'infinity'  downstream  in  the  wake).  This  Joe a have  a marked  'bucket', 
due  to  the  rapid  movement  of  transition  on  the  upper  surface,  back  from 
x/c  " 0.09  at  M^  «■  0.72  (Fig  11)  to  O.b  at  - 0.745  (Fig  7).  Without  this 
transition  movement,  the  value  of  the  total  drag  coefficient  calculated  at 
Mv  - 0.745  would  have  been  much  higher,  by  about  0.004;  so  that  drag-rise  Mach 
number,  , would  have  been  predicted  as  about  0.73,  as  opposed  to  the  apparent 

value  of  0.75  from  Fig  lb. 

A pictorial  summary  of  the  comparisons  of  pressure  (Mach  number ) distribu- 
tion is  provided  by  Fig  17,  which  gives  on  a small  scale  a 'traverse'  through 
the  Mach  number  range  O.b  to  0.77  at  constant  “ 0.85°,  and  through  the  range 
of  a from  -2°  to  *2l  at  M “ 0.745. 

IV 

5 PAUHILATIONS  AT  H1CHKR  REYNOLDS  NHMBKKS 

Because  of  certain  unsatisfactory  features,  referred  to  above,  of  the 
experimental  information  at  present^  available  for  this  aerofoil  - unavoidable 
at  the  low  Reynolds  number  of  the  Nl.R  tests  - it  was  decided  to  make  a purely 

* See  footnote  on  page  5 

**  In  the  sense  that  they  would  not  occur  in  fully  turbulent  flow,  even  at  high 
Reynolds  numbers  (see  section  5). 

t The  aerofoil  has  now  been  tested  in  the  Lockheed  Georgia  Compressible  Flow 
facility  at  high  Reynolds  numbers,  but  results  are  not  yet  available. 


t heore  t i i*4 1 study  of  its  behaviou:  .it  two  highci  Reynolds  numbers:  10  , being 
tvpical  o!  what  would  be  obt  ainab  1 e in  a laiger  wind  t unne  1 , and  S ' 10’,  a 
typical  full-scale  value.  The  present  method  was  used,  with  transition  fixed  at 
0.0'  chord  on  both  surfaces,  first  'traverses’  in  Mach  number  were  made  at  con- 

7 

stant  values  of  0 - 0.  ',  O.s,  0.S,  0.SS  iRe  - 10  only),  O.b  and  0.b5 

7 ^ 

tRe  » S " 10  onlv)  - starting  in  each  case  at  M^  - O.b  and  increasing  in 

sufficiently  small  steps  to  enable  the  drag  rise  Mach  number  iM^)  and  the  Mach 

numbers  tor  onset  ot  trailing  edge  separation  tM  ) to  be  determined  with  reason* 

sop 

able  accuracy.  The  former  iH^)  is  here  convent i onal ly  defined  to  be  the  value  ot 

M tor  which  0 is  0.00.0  greater  than  its  value  at  M * O.b  at  the  same  lift 

coefficient:  except  that  tot  0,  ' 0.S  there  is  some  wave  drag  i ;’<•  C 

l D total 

i so  that  C„  at  M - 0.6  has  then  been  used  as  the  datum  value, 

viscous  P viscous 

The  latter  tM  ) has  been  estimated  bv  using  the  empirical  observation  that 
s e p 

appreciable  ettects  ot  trailing  edge  separation  fin  particular,  divergence  of 
trailing  edge  pressure'  are  found  in  practice  when  the  transformed  boundary  layer 
shape  parameter  H exceeds  about  2.2*  one  boundary  layer  thickness  upstream  of  the 
trailing  edge  on  the  upper  surface,  U1  ^ ^ ' 2.2).  Next,  in  order  to  extend  the 
estimated  ’separation  boundaries’  to  lower  Mach  numbers  and  higher  lift  coef- 
ficients, traverses  in  a at  constant  values  ot  M were  made,  increasing  a 
until  separation  conditions  were  obtained,  as  defined  above.  It  was  found  however 
that  this  procedure  was  only  satisfactory  at  Mach  numbers  of  0.72  and  0.70;  for 
lower  values  ot  M the  shock  waves  became  so  strong  that  excessive  values  o t 


h r 


2.S)  were  predicted  at  the  foot  of  the  shock  well  before  H ^ reached  the 


critical  value.  In  such  c i re  unis t ances  the  prediction  of  separation  must  be  treated 
with  reservation  in  our  present  state  ot  knowledge. 

before  examining  the  predictions  of  M and  M , we  show  some  typical 

p sop 

pressure  distributions.  First,  for  He  - 10',  Fig  ISa.b&o  show  the  results  for 

i',  “ d . ■< , for  M from  O.b  through  0.  ’•«  fM  ) up  to  0.7b  tbevond  M '.  Similar 
1 •'  P sep 

curves  for  0,  ■ 0.S  are  shown  in  Fig  l'>a.tb  iM,  ' 0.7'.  M ' 0.7s',  and  for 
1 ' 0 sep 

0 ■ O.b  in  Fig  20a&h  ' O.t'S,  M ' 0.7,').  Corresponding  drag  curves,  tor 

i.  i»  sop 

Oj  “ O.s,  0.S  and  O.b,  are  shown  later  in  Fig  2b.  Fig  21  shows  how  the 
pressure  distribution  varies  with  0,  for  M - O.bx;  at  the  highest  value  ot 
Cj  , '.02,  separ a t i on  is  predicted  at  the  trailing  edge  but  would  undoubtedly 
have  occurred  earlier  at  the  foot  ot  the  shock.  Temporarily,  this  situation  can 
lead  in  practice  t o an  appreciable  m l i 1 1 coefficient,  as  long  as  the 
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* It  is  important  t o note  that  tltis  particulai  value  should  onlv  be  associated 
with  the  present  method:  it  is  not  intended  t o implv  anv  absolute  significance 
to  it. 
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separation  bubble  does  not  reach  the  trailing  edge;  so  that  the  predicted  value 

tor  Ci  mav  not  in  fact  be  optimistic. 

‘•sep 

Turning  now  to  the  calculations  at  full-scale  Reynolds  number  (5  » 107), 

we  show  in  Kig  22a&b  how  the  pressure  distribution  varies  with  M at  C “ O.b; 

and  in  Kig  23aSb  at  C - 0.b5.  At  first  sight  the  effect  of  the  increase  in 

7 

Reynolds  number  appears  small;  but  a close  comparison  of  Kigs  20b  [Re  • 10  , 

C|  ■ O.bl  and  22b  [Re  “S'  107,  C.j  “ O.bl  shows  that  the  shock  strength  is 
slightly  weaker  at  the  higher  Reynolds  number.  Some  pressure  distributions  at 
- 0.72  with  Cj  increasing  up  to  separation  conditions  (C^  - 0.83)  are  shown 
in  Kig  24. 

An  interesting  comparison  is  made  in  Kig  25,  in  which  are  shown  at  the 

'design  point'  IN  * 0.72,  C * 0.6)  the  pressure  distributions  at  Re  “ 10*  and 
7 * L 

5 > 10  , compared  with  the  shock-free  distribution  given  by  the  full  line.  At 
Re  - 10  a shock  is  present  sufficiently  strong  to  cause  appreciable  wave  drag 
iso  that  this  point  is  well  beyond  the  drag-rise  boundary  - see  Kig  27);  but  at 
Re  » 5 *<  107  the  shock  is  weaker  and  in  fact  this  condition  turns  out  to  be  an 
optimum  one  (giving  the  maximum  value  of  Mh/D)  at  full  scale  for  this  aerofoil. 

Some  typical  curves  showing  the  variation  of  C with  Mw  at  constant 
are  shown  in  Kig  2b.  At  • 0.4  there  is  little  scale  effect  on  either  drag 

creep  or  . At  0^  » 0.5  there  is  much  less  drag  creep  full-scale  than  at 

Re  ■ !07,  though  is  not  greatly  affected.  But  at  * O.b  the  drag  creep 

is  so  large  at  Re  * 107  that  is  only  0.64,  as  opposed  to  0.725  at  Re  * 5 ' 10*. 

The  final  figure.  Fig  27,  gives  the  drag  rise  and  separation  boundaries 
estimated  for  this  aerofoil  at  the  two  Reynolds  numbers.  As  mentioned  above,  a 
strong  favourable  scale  effect  on  is  predicted  for  > 0.5;  and  as  regards 

the  separation  boundary  the  effect  is  favourable  at  all  lift  coefficients  (above 
0.4),  though  it  does  increase  further  with  . Thus  under  full-scale  conditions, 

the  aerofoil  at  its  optimum  condition  « 0.72,  Cj  * O.b  should  have  a margin  to 
separation  of  about  0.02  in  Mach  number  and  0.23  in  lift  coefficient. 

We  have  also  included  in  Kig  27  sketches  of  the  Mach  number  distribution  at 

a number  of  points  on  the  drag-rise  and  separation  boundaries  (full-scale), 

showing  how  the  shape  of  the  curves  varies  over  a wide  range  of  conditions.  It 

is  interesting  to  note  that,  at  the  optimum  point  (M  - 0.72,  C,  « O.bl,  the 

L 10 

pressure  distribution  corresponds  closely  to  Whitcomb  s reconmendat ions  for  a 
good  supercritical  aerofoil,  with  a weak  expansion  just  behind  the  shock  wave: 

119  but  with  the  added  advantage  of  a more  'peaky'  upper  surface  pressure  distribution. 
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with  appreciable  isentropic  compression  (from  M—  l.J  to  about  1.2)  ahead  of  the 
shock . 

Thus  the  aerofoil  shows  every  sign  of  being  an  extremely  good  one  in 
practice;  but  as  we  have  mentioned  earlier  the  existing  experimental  information 
about  it  is  extremely  misleading.  For  example  Fig  lb  would  indicate  a value  of 
My  of  about  0.755  (at  ■ 0,i5);  this  point  has  been  marked  on  Fig  27,  and  it 
can  be  seen  that  this  value  of  is  about  0.02  higher  than  the  value  predicted 

fUli-saalti  for  the  same  . It  must  be  emphasised  that  this  does  not  appear 

to  be  a blockage  problem:  it  is  entirely  due  to  the  suppression  of  the  shock 
wave  at  the  low  Reynolds  number  of  the  tests  by  interaction  with  the  laminar 
boundary  layer. 


1.1  ST  OF  SYMBOLS 


aerofoil  chord 
drag  coefficient 

theoretical  viscous  drag  coefficient  (.from  wake) 

lift  coefficient 

pressure  coefficient 

boundary  lover  shape  parameter 

free  stream  total  pressure 

Mach  number 

free  stream  Mach  number 
drag  rise  Mach  number 
static  pressure 

Reynolds  number  (based  on  chord) 
velocity  at  edge  of  boundary  layer 
('artesian  coordinates  (see  Fig  11 
value  of  x c at  transition  point 
ang loot  i n>  i dence 

angle  of  incidence  in  wind  tunnel  (uncorrected) 

irtificial  viscosity  parameter 

shock  parameter 

kinematic  viscosity 

boundary  layer  momentum  thickness 
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